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Abstract — The mechanisms for the formation of X-ray lines in the spectrum of SS 433 are investigated by taking 
into account the radiative transfer inside the jets. The results of Monte Carlo numerical simulations are presented. 
The effect of a decrease in line intensity due to scattering inside the jet turns out to be pronounced, but it does not 
exceed 60 % in magnitude on the entire grid of parameters. The line broadening due to scattering, nutational motion, 
and the contribution of satellites can lead to overestimates of the jet opening angle G from the line widths in Chandra 
X-ray observations. The fine structure of the lines turns out to be very sensitive to the scattering effects. This makes 
its investigation by planned X-ray observatories equipped with high-resolution spectrometers (primarily Astro-H) a 
powerful tool for diagnosing the parameters of the jets in SS 433. 
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1 INTRODUCTION 

The Galactic source SS 433 is currently the only close 
X-ray binary system observed in a permanently supercriti- 
cal regime of accretion. This manifests itself most clearly 
in the presence of a pair of oppositely directed relativistic 
jets oriented perpendicular to the plane of a s upercritical 
accretion disk (for a review, see Fabrid (l2004b ). 

X-ray observations of the system show that apart 
from the continuum component attributable mainly to 
bremsstrahlung from the jets, its spectrum exhibits numer- 
ous lines of highly ionized atoms of heavy elements. By 
measuring their characteristics, one can get an idea of the 
physical conditions in the regions closest to the compact 
object in which the jets ar e formed, accelerated, and colli- 
mated (IKotani 



:n tne jets ar e 
i et al.lll996 l). 



Substantial progress in understanding these fun- 
damental mechanisms has been achieved in the 
last ten years throu gh observation s by the Chandra 



INTEGRAL( Cherepashchuketal 



2005h. 



dMarshall et al.l 120021: Ikamiki et^ 120031: iLopezetaL. 
l2006h . XMM-Newton dBrmkmann et al.1 l2005h . 



RXTE 



dFilippova et al. 200i), and Suzaku dKubota et al.ll201Qh 



X-ray observatories in the standard and hard spectral 
ranges. However, the uncertainty in physical parameters 
of the jets (density, size, and opening angle) remains 
significant (for more details, see Sections 4.3 and 6). With 
high- spectral resolution instruments (primarily Chandra 
HETGS), it has been possible to measure the line widths 
and intensities and to estimate the opening angle, density, 
and temperature of the hottest parts o f the j et s from them 
dMarshall et"aD (|2002),Namiki etaD d2003L iLopez ( ' 



: et al.1 



2006h ). However, these required correcting the previously 
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adopted standard jet model (IKotani et al. Il996h . For 
example, the emission measure of the hottest parts of the 
jet estimated from the line intensity in the high-energy part 
of the spectrum turned out to be approximately half the 
value expected from the standard model. In addition, it was 
required to assume an overabundance (about 30% relative 
to the solar value) of iron, sulfur, silicon, magnesium, and 
neon ([Marshall et al. 1 12002 ) for the observed continuum 
to agree with that predicted from the line intensities. The 
EPIC/XMM-Newton spectra point to a more than eightfold 
excess of radiation near the Kg triplet of helium-lik e nickel 
dBrinkmann et al.l d2005l)rMedvedev. Fabrikal d2QTa) ). The 
physical interpretation of such discrepancies was often 
reduced basically either to the uniqueness of the SS 433 
phenomenon or to the inclusion of additional emitting 
or absorbing components of the system while restricting 
oneself to using a very simple scenario for the formation 
of X-ray lines. 

In this paper, we study in detail the influence of scat- 
tering effects inside the jets on the spectrum of emergent 
radiation in lines of the standard X-ray band for the first 
time. The necessity of such a study follows from simple 
estimates of the jet transverse optical depth for Thomson 
scattering (~ 0.1) and, what is especially important, for 
resonant scattering on electron transitions in highly ionized 
atoms of heavy elements (~ 100 for the permitted tran- 
sitions in helium-like ions). As a consequence, the effect 
of an increase in the mean free path of a photon before its 
escape from the jet due to multiple resonant scatterings be- 
comes important. This, in turn, leads to an increase in the 
probability of being scattered by a "hot" electron and leav- 
ing the line while reducing its intensity and forming broad 
wings. 

Using the Monte Carlo technique for radiative transfer 
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problems developed by IPozdnyakov et aP (1983) allowed 
the formation of X-ray lines in the spectrum of the jets to 
be simulated most accurately by taking into account the 
multiplicity of the most important of them and the possi- 
ble deviations from the coronal approximation. As a result, 
we made predictions regarding the observational line char- 
acteristics, bearing in mind the high resolution of planned 
spectroscopic experiments (primarily Astro-H; see, e.g., 
Irakahashi et al.,(.201Q) . 



2 THE JET MODEL 
The standard multitemperature jet model for SS 433 



(iKotani et al. 1996; Mars hall et a l. 2002) served as a start- 
ing point for our study. In this model, the X-ray jet is 
treated as a highly collimated (an opening angle 6 1° 
) ballistic plasma flow directed away from the compact ob- 
ject perpendicular to the plane of a supercritical accretion 
disk. The velocity component parallel to the jet axis is con- 
stant at each point and equal to 0.26c. The matter is dis- 
tributed uniformly within each layer perpendicular to the 
jet axis. Under the assumption of an axisymmetric flow 
pattern, it is convenient to introduce a conical coordinate 
system^ (r,6>), where r is the distance measured from the 
cone vertex along the jet axis and 6 is the angular displace- 
ment from the axis. The directly observable jet region clos- 
est to the compact object will be called the jet base. The 
location ro, electron density rieo, and temperature To of this 
region as well as the opening angle gB are the input param- 
eters of the model. 

Given the uniformity of the matter distribution in the 
r = const layer, the law of change of the electron density 
follows from the condition for the conservation of mass 
flux along the jet: 



(1) 



Jo 



Using the derived density profile, we can write the ther- 
mal balance equation as 



dr r 



lUeUi A(T) 
3(ne-\-ni) V 



(2) 



where the first and second terms on the righthand 
side correspond to the cooling through adiabatic expan- 
sion and the losses through radiation, respectively. Intro- 
ducing dimensionless quantities 7^=T/To, ^=r/ro, X^n/Zn^, 
A23 (^)=A (vTo)/(10-^^^^^) and assuming that 7=5/3, we 
obtain 



dr]___Ar]__ A23 ij]) 



(3) 



* In what follows, the temperature is in energy units and the opening 
angle is in radians. In this case, since the latter is small, we assume that 
tan Q ^ Q. 



a =10 



-23 



erg cm-^ 2 Ueoro X 



(4) 



s 3 vTo 1+X 

Substituting the input parameters typical of the jet in SS 
433 and assuming that X ^ 0.91, we find 



a ^0.1272 



Vl0i4cm-V VlOiicm/ V20keVy V0.26cy 



(5) 



The limits of small and large a correspond to adiabatic 
and quasi-cylindrical jets, respectively (see below and Sec- 
tion 4.3). 

The emissivity A(T) in Eq. ^ is calculated in the 
model of a ho t, optically thin plas ma using the low-density 
limit (APEC, ISmith et al.' ("200 1")) by assuming solar ele- 
mental abundances ((Lodders 2003)). A numerical solu- 
tion of Eq. ([3]) with the parameter a and the initial con- 
dition r]a 1^=1= 1 gives the temperature profile along the 
jet r]a (0 or T(r) = Tor]a (^/^o)- Since r]a (0 is a mono- 
tonic function of ^, there exists an inverse function ^q, (t]). 
This allows the boundary of the computational domain ^^ax 
to be determined, because thermal instabilitie s emerge in 
the jet at Tmin ^ 0.1 keV dKotani et al.lll996h : £,max(o^) = 

The number of free parameters can be reduced by fixing 
the total X-ray luminosity: 



where 



I(a) = 7rX 



I 



d^ 



(6) 



(7) 



The cooling at low a is determined by adiabatic expan- 
sion, i.e., r]a(0 ^ ^o(0 = r"^^^ and £,max = (Tmin/To) ^^"^ 
do not depend on a. Consequently, 1(a) ^ const and 
Lx = const implies 

n^Q rl = const . (8) 

At large a, in view of Eq. © AiZMi:^ ^ ^ Con- 
sequently, 1(a) (xl/a and = const implies 



rieo ^0 



rn 6^ = const, 



(9) 



which also corresponds to a constant mass loss rate in the 
jet My = const. 

3 RADIATIVE TRANSFER IN LINES 

It is convenient to investigate the radiative transfer in a 
comoving frame of reference by getting rid of the constant 
longitudinal velocity component. The velocity field in such 
a frame is 

V|| (^) =0, vi_ (^) =Vj^tan(0), = 0.26c (10) 
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In the introduced frame of reference, we will consider 

a line photon with energy E emitted on the jet axis in a 
— ^ 

direction perpendicular to the jet axis. Since E^T <^ 
nieC^ = 511 keV in the situation of interest to us, the optical 
depth to the jet edge for scattering by free electrons for 
such a photon is rr(0 = ^eCO^r^o©^, where cft = 6.65 • 
10~^^ cm^ is the Thomson cross section. The calculation 
of the optical depth to the jet edge for resonant scattering 
on electron transitions in ions is slightly complicated by 
the velocity field (fTOb . Therefore, we will initially assume 
that E = Eq, where Eq is the energy of the corresponding 

transition. Given the Doppler shift in the velocity field ([TOb 

— y 

and the fact that the direction Qq at each point is along the 
local velocity of the matter, we will then obtain 



v±(0)Eo Y 
c AEdJ 



r^idO, (11) 



where G [0, 6] and, in view of the smallness of 6, along 
with 0, we set tan (6)^6. We also used the following stan- 
dard notation: ctq is the resonant scattering cross section at 
the line center and AEd is the line width due to the thermal 
motion of ions (for more details, see Section IT. 1.21 ). 

Denoting x=f^(l3 = vj/c = 0.26), y = x^, Tres (0 = 
(0 cro^o©^ and integrating over we will obtain the total 
optical depth to the jet edge for resonant scattering for a 
photon emitted on the axis in the direction Qq: 

te-y'dy 

^(x) 

^effiO = ^resiO = ^re^(0 , (12) 

XX 

where ^(x) = ^erf(x), erf(x) is the standard error func- 
tion. 

Let us now consider a photon also emitted on the jet axis 
in the direction Qq but with an initial energy shift A£^o rel- 
ative to the central value £"0, i.e., £ = £"0 + AEq. We see that 
for such a photon 

dTres(0,0 = ni(OcToey ^ r^^dO , (13) 

and the total optical depth to the jet edge 



/ e-y dy 



Teff(OAEo) = rres(0- 



(14) 



■ rresiO 



^((l-e)x) + ^(ex) 



where e = is the initial shift parameter (in addition, 
we used the oddness of erf(x)). Thus, the effective optical 
depth for resonant scattering can generally be written using 
the introduced notation as 



The function tz, (x) = (see Fig. [B is such that: 

1) K,e(x) < I for any x and e; 

2) K,e(x) is symmetric relative to e = 0.5, i.e. K,o.5-5e(x) = 
1^0.5+6 e (x) for any x and ^e; 

3) for any x and e 72^^0.5 hz^ (x) < /^o.5 (x); 

4) at large x, hzo (x) ^ l/x; 




Fig. 1. General view of the function ^Ce(x), characterizing the resonant 
scattering efficiency as a function of x = and the shift parameter 

e = according to Eq. (T5\ 



Thus, it can be concluded that not the entire optical depth 
Tres but Only some part of it Te/f is involved in the resonant 
scattering of the photon under consideration; the smaller 
the ratio of the thermal ion velocity to the bulk one speci- 
fied by the velocity field (fTOl) , the smaller this part. In view 
of the asymptotics defined by property 4), this effect for a 
photon at the line center (e = 0) is interpreted in such a way 
that the scattering at large x occurs only inside the cylinder 
at the boundary of which the bulk ion velocity is equal to 
the thermal one. This means that the estimate of re//(^,0) 
remains valid not only for the photons emitted on the jet 
axis but also at any other point, except for the surface layer 
of thickness roQ^/x. In addition, in view of properties (2) 
and (3), we can assume the profiles to be asymmetric rel- 
ative to the central value due to the larger effective optical 
depth for high-energy photons (E > Eq^c > 0) (as an illus- 
tration of these properties, see Fig. O. 

Let now the distribution of sources along the jet be de- 
fined by some function ^(0^ i-^-^ the relative contribution 
of the layer d(_ to the total photon luminosity in the line is 

(j)(OdO and J ^(0<^^ = 1- The weighted mean optical 



depths 



(15) 



Tt 



Cmax 



(16) 
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Fig. 2. Resonant scattering of photons in a layer perpendicular to the jet 
axis for x = = ^-0- ~ 5. 1,2,3 - the pattern of the trajectory for 

photons emitted on the axis with E = Eq^E > Eq^E < Eq respectively; 4 
- the pattern of the trajectory for a photon emitted not on the axis with 
E = Eo 



4.1.1 Emission 

We will be interested in the radiation from a jet in the line 
corresponding to the transition from the upper level (UL) to 
the lower level (LL) in an n-fold ionized atom of element 
Z (for example, for the resonance Kq, lines of heliumlike 
atoms n = Z-2, UL= ls2p(iPi) and LL = Is^C^^o). 

The dependences ne(^) and T found above allow 
the local ionization balance to be computed in the coronal 
approximation using AtomDB 2.0. As a result, in each 
^ = const layer of thickness d^ (J^ <C 60 we will obtain 
the density of ions of a given type ^/(O ^nd the plasma 
emissivity JiT) in the line of interest to us at temperature 
T = Tor]a (0- The line photon production rate in such a 
layer is then defined as 



dL(0 = 7rnl,rle^ 



e 



Thus, <l^(0 = T^^^ where the total luminosity 



(18) 



'^eff 



^max 

-I 



(17) 



then characterize the influence of scattering on the inte- 
grated line emission. 

To estimate the probability of the scattering by a free 
electron, we can us e the critical op t ical d epth r^r = 
following IPozdnyakov et al.l (Il983h . At tt < 



2 Inijeff/TT) 

Tcr < 1, the jet may be considered optically thin for scatter- 
ing by electrons even for resonant photons. At r^r < < 1, 
almost all of the photons will be scattered by the electron 
only once and, having left the line, form broad wings. 

Undoubtedly, these estimates are qualitative in nature 
due to the dependence of Te/f on photon energy, the distri- 
bution of sources inside the layer d^, the velocity field ([TOb , 
and the line overlap in the case of doublets and triplets. For 
a detailed analysis of the line formation, we have to resort 
to numerical simulations. 

Since this work is aimed not at a theoretical study of 
the radiative transfer in lines but rather at its application 
for a specific physical object, we will restrict ourselves to 
the presented analysis, because the estimates obtained in 
its context can give a clear interpretation of the results ob- 
tained in our simulations (see the "Results" Section). 



4 SPECTRUM SIMULATION 
4. 1 The Scheme of Computation 
The approach used was developed by lPozdnyakov et aD 



( 19831) and, subsequently, was successfully applied 
to investigate the radiative transfer in g alaxy clusters 
( Sazonov et aDl20Q2l: Izhuravleva et al.ll2QTQl) . 



Lo = 7rnl,rle' J 



1 



^2 ^S- 



(19) 



Similarly, for the continuum radiation 



Lc(E) = 7rn,Qr^e ' 



I 



e 



di, 



(20) 



where JdE^T) is the plasma emissivity in the contin- 
uum per unit energy interval with center E at temperature 
T computed using the NoLine model of AtomDB 2.0.1 
(kindly provided by Adam Foster [3). The equivalent width 
of a line with a transition energy Eq is then defined as 



EW = Lo/Lc(Eo). 



(21) 



However, it is worth noting that the equivalent width 
defined in this way has no direct bearing on the observed 
quantities, because both jets contribute to the actual contin- 
uum, given the Doppler shift and relativistic collimation of 
the radiation. Nevertheless, assuming the jets to be iden- 
tical, determining the observed equivalent width from the 
line equivalent width that we use for an arbitrary precession 
phase does not seem problematic. 

4.7.2 Scattering 

The reciprocal of the mean free path (scattering coeffi- 
cient) serves as a local characteristic defining the photon 
scattering probability: 



Er = nrO-r 



(22) 



* http://atomdb.org/; 

"^A description of the analogous model for APEC v 1.3.1 can be found 
at |http://cxc.harvard.edu/twiki/bin/view/SnrE0102/NoLine| 
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where % is the number density of scattering centers and 
(Jx is the cross section for the corresponding scattering. For 
example, for scattering by electrons, = UeCFT, cft = 6.65 • 
10~^^ cm^ is the Thomson cross section. For resonant scat- 

2' 



tering by ions = riiares, Ores = ctq exp 



2T 



Am„ 



1/2 



\^Ed ) 



l\E=E-EQ{l + ^-fcos^lj), 
where is the classical electron radius, nip is the proton 
mass, A is the atomic weight of the ion, E is the photon en- 
ergy, and \Ij is the angle between the photon direction and 
the local gas velocity determined by field ([TOl) . The energy 
£"0 and oscillator strength / of the atomic transition can be 
found using AtomDB 2.0.1. It is worth noting that although 
the oscillator strength / is a basically positive quantity, in 
what follows (e.g., in Table 1), we take / = for transitions 
with / < 0.0001 to simplify our calculations. 

Whereas depends only on the local electron density, 
is a function of the local plasma characteristics {rie and 
r), the photon energy and direction. 

We performed detailed Monte Carlo simulations of the 
emergent radiation spectrum in lines. We used a scheme 
with a statistical "weighting" of photon pack ets, along wit h 
the method of a constant total cross section (ISoborlll973l) . 
The idea of this method is to introduce, along with the real 
types of scattering, some fictitious scattering, such that 



^fict = So - (Se + S/) , 



(23) 



where the constantSo is greater than or equal to the max- 
imum possible T^reai = Sg + E/. As a result, the inhomo- 
geneous (in the sense of local scattering characteristics) 
medium is replaced by a homogeneous one with a small 
mean free path A = In each scattering, the realization 
probability of a particular type is defined as px = ^x/^o- In 
the case of fictitious scattering, both direction and energy of 
the photon package remain unchanged. The resonant scat- 
tering by ions is represented as a combination of the dipole 
and isotropic components. The weight of the dipole com- 
ponent Wd is determined by the total angular momentum 
of the lower leve l Jll and by the difference Aj = Jul - Jll 
(lHamiltonlll947l) . 

The described scheme can be extended to the simulation 
of unresolvable multiplets (the photons fall into the total 
spectrum "weighted" proportionally to the parent line in- 
tensity, while the scattering on each multiplet component 
is considered as a separate type of scattering). The addition 
of a broad continuum component and the calculation of its 
scattering are also possible. 

4.2 The Set of Lines 

To form the set of simulated lines, the best res- 



olution Chandra HETGS spectra (iMarshall et al.l 12002 : 



iNamiki et al1l2003l: iLopez et al.ll2006h should be used. The 
brightest observed lines subject to strong scattering effects 
correspond to the transitions in helium-like ions and to 



the Lya transitions in hydrogen-like ions with significant 
oscillator strengths whose parameters are listed in Table [T] 
In this case, the transitions in helium-like ions have a 
triplet structure (the resonance (w), intercombination (x + 
y), and forbidden (z) lines). Bearing in mind the impor- 
tance of these lines for diagnosing the jet plasma parame- 
ters, we performed detailed simulations of the components 
by taking into account the possible interaction via the opti- 
cal depths and the influence of collisions on the upper level 
population of the forbidden and intercombination transi- 
tions at ne greater than some ng^crit (Porquet et al. 2oTo|). 
In addition, in some cases, the influence of satellites on the 
intensity ratio of these lines should be taken into account 
(see the "Results" Section). 

4.3 Input parameters 

There is a considerable uncertainty in the physical pa- 
rameters of the jets in SS 433. Based on Chandra HETGS 
observations at the phase of the gr e atest d isk opening to- 
ward the observer, Marshall et al.l (l2002b obtained ro ~ 



2 X 10^0 cm, neo^2x 10^^ cmT^ Tq ^ 13 keV, and 6 ^ 
0.01 rad for the approaching jet. Analysis of the XMM- 
Newton X-ray spectra also at the phase of the greatest 
opening dMedvedev. Fabrikall2010h gives ro — 2 X 10 cm 



and To 17 keV for the approaching jet. In contrast, stud- 
ies based on the eclipse of the jets by the op tical companion 
yield ro ~ 1 x 10^^ cm and Tq ~ 30 keV dFilippova et al. 
2006h . 

Therefore, two classes of models satisfying condi- 
tion © for the conservation of total X-ray luminos- 
ity were formed - for t he pair of jets Lx ^ 10^^ erg /s 
(iMedvedev. Fabrikall2010l) i n the ranse from 0.1 to 50 keV 
either in form ^ (the class of quasi-adiabatic models) or in 
form © (the class of quasi-cylindrical models), depending 
on the parameter a. Since the available observations can 
reliably set only an upper limit on the opening angle (see 
Section 6), we allowed 6 to change in a very wide range. 
For this reason, precisely 6 served as a class-forming pa- 
rameter (in the sense of discrimination by a). As a result, 
the simulations were performed on a wide grid of input 
parameters given in Table |2] and encompassing the most 
probable values. The quasiadiabatic models correspond to 
6 = 0.01 rad and O = 0.02 rad, while the quasi-cylindrical 
ones correspond to 9 = 0.003 rad and 6 = 0.0001 rad. We 
considered ro=20 keV and 7o=30 keV as possible plasma 
temperatures at the jet base. The abundances of heavy ele- 
ments were assumed to be solar (Lodders 2003) . 

It is worth noting that whereas in the case of quasi- 
adiabatic models ro roughly coincides with the distance 
from the jet base to the compact object, in the case of quasi- 
cylindrical models ro has the meaning of only a geometri- 
cal parameter of the flow. Therefore, using the combina- 
tion ro6, which means the transverse size of the jet near 
its base, appears more preferable for the description of a 
particular model in the general case. 
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Table 1. The set of simulated lines. The asterisk (*) marks the lines for which we do not present our simulation results but which are used in constructing 
the total broadband spectrum (see the Conclusions). The expression / = should be understood in the sense that / < 0.0001 (see also the text) 



Spectroscopic 
syniDOl 


Transition 


Eo, keV 


f 




FeXXV 


: \s^-\s?>pCPi) 


7.881 


0.14 


1 


Ni XXVII 


Ka : ls^-ls2pCPi)(w) 


7.806 


0.72 


1 


Ni XXVII 


Ka : ls^-ls2p CP2)(x) 


7.799 





- 


Ni XXVII 


Kc,'As^-ls2p (M)(y) 


7.766 


0.07 


1 


iNl AAV 11 


Ka . is - VSlp ( Oij(Zj 


n HAA 
1 . /44 


U 




Fe XXVI 


Lj,: 1^-2/7 (1^3/2) 


6.973 


0.27 


0.5 


Fe XXVI 


Lyo,-As-2pCP,,2) 


6.952 


0.14 





Fe XXV 


Ka : ls^-ls2pCPi)(w) 


6.700 


0.78 


1 


Fe XXV 


Ka : ls^-ls2p CP2)(x) 


6.682 







Fe XXV 


K^'As^-ls2p (M)(y) 


6.667 


0.07 


1 


Fe XXV 


Ka : ls^-ls2p CSi)(z) 


6.636 







Ca XIX^ 


Ka : ls^-ls2p CPi) 


3.902 


0.77 


1 


Ar XVIl3 


Ko^'As^-\s2pCPi) 


3.133 


0.77 


1 


SXVI^ 


Ly^As-2pePy2) 


2.623 


0.27 


0.5 


sxvi3 


Ly^As-2p(}P,,2) 


2.620 


0.14 





SXV 


Ka : ls^-ls2p ^Pi)(w) 


2.461 


0.76 


1 


SXV 


Ka : ls^-ls2p CP2)(x) 


2.449 







SXV 


KaAs^-ls2p (M)(y) 


2.447 


0.07 


1 


SXV 


Ka : ls^-ls2p CSi)(z) 


2.430 







Si XIV ^ 


Lyo,As-2pePy2) 


2.006 


0.27 


0.5 


Si XIV 3 


Lyo^As-2pCPi/2) 


2.004 


0.14 





Si XIII 


Ko^As^-\s2pCPi){w) 


1.865 


0.75 


1 


Si XIII 


Ka : ls^-ls2p CP2)(x) 


1.855 







Si XIII 


K^As^-ls2p (M)(y) 


1.854 


0.07 


1 


Si XIII 


: ls^-ls2p CSi)(z) 


1.839 
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To illustrate the physical picture that emerges when 
considering the constructed grid of parameters through the 
prism of the jet model from Section 2, it is convenient to 
use the profiles of temperature T (0 (Fig. O, optical depth 
(0 = e^roS;,, and flux (j)i (0 (Fig. H in lines. 









= 0.01 rad 






To = 20 keV 




















1 ?° \ 


\ CP \^ 






1 ^ \ 






1 


\ % 

\ <^ 

a 






Vj/c = 0.26). The results for the remaining (marked by * 
in Table [B lines are used only in constructing the total 
broadband spectrum (see the Conclusions). Analysis of our 
computations for other precession phases and To = 30 keV 
showed that the general conclusions (see the Conclusions) 
remain valid here as well, while most of the changes are 
equivalent to some shift on the constructed grid of parame- 
ters in the sense of scattered line characteristics. 

The main quantitative characteristic reflecting the influ- 
ence of scattering effects on the observed spectrum is the 
fraction of the photons C that left the line as a result of their 
scattering by a free electron and that fell into the broad 
wings. Together with the temperature Te providing the 
best agreement between the shape of th e broad wings and 
the s ingle Compton scattering kernel (ISazonov. Sunyaev 



2000), the parameter ^ allows the contribution of scattered 



Fig. 3. Temperature profiles along the jet for various input parameters 
from TableE] 



radiation to the continuum to be completely described (for 
more details, see Section 5.2). Since it is probably impos- 
sible to find such a complete and universal way of describ- 
ing the shape of the scattered line itself, one has to content 
oneself with standard approximations when analyzing the 
derived profiles. This approach is justified from the view- 
point of observational characteristics of existing spectrom- 
eters (Chandra HETGS), which do not allow the fine struc- 
ture of the line profiles to be investigated in most cases. 
However, it is clearly insufficient from the viewpoint of 
new-generation X-ray observatories (primarily Astro-H). 
For this reason, the set of parameters, except for the line 
attenuation coefficient equal to (1 - Q, is determined by the 
specific simplifying model that is most appropriate for the 
line in the resolution of a particular instrument (Chandra 
HETGS (to be more precise, HEG-High Energy Gratings) 
or Astro-H SXS). 

To analyze the spectra in the Chandra HEG resolution, 
we used a model response function in the form of a Gaus- 
sian with FWHMchiE,) = 29.9 eV {^j^f (in accor- 
dance with the Chandra Gratings Fact Sheet Q), where 

= Ez/(l +z). For Astro-H SXS, we used a Gaussian 
with FWHMah (Ez) = ^eV = const (in accordance with the 
Astro-H Quick Reference Q). 

Since all lines exhibit very similar patterns in various 
regimes of radiative transfer, we will restrict ourselves to 
a detailed description and interpretation of the results only 
for the Ka triplet of helium-like iron, providing the scatter- 
ing features in a specific line, where necessary. 



5 RESULTS 
5. 1 Lines 

To save space, we present the results of our compu- 
tations only for the most important (in our view) lines of 
the approaching jet with To = 20 keV at the phase of the 
greatest disk opening that corresponds to an angle between 
the jet axis and the observer's direction x ^ 65° , i.e., a red- 
shift =-0.078 (z= 7(1- /^^^^X)-l, 7=l/v/l-/^^ f^ = 



5.1.1 FeXXVKa 

The helium-like iron triplet is brightest among the ob- 
served lines (Table O. Nonetheless, only the components 
corresponding to the allowed transitions with / > (see 
Table [B can have a significant intrinsic optical depth for 

* https://icxc.harvard.edu/rws /peer_review/proc_docs/ 
/Fact_Sheet_gratings .pdf 

"^http://astro-h.isas.jaxa.jp/doc/ahqr.pdf 
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resonant scattering. At the same time, the Hne broadening 
due to the transverse velocity component at 6 ~ 0.01 rad 
is enough to overlap the resonance and intercombination 
lines. This makes the interaction between the components 
via the optical depths possible and, consequently, indicates 
that all components should be taken into account simul- 
taneously when the local characteristics of the scattering 
medium are calculated. 

The computation perf ormed with th e Clou dy code (ver- 
sion 08.00 described bv iFerland et al.l (Il998)) in the tem- 
perature range of interest to us showed that for helium-like 



iron Uecrit ~ 10^^ cm ^ is the density at which the Colu- 
sa 



sional excitations from ^Si (the upper level of the forbid- 
den component) to ^^^3,1,2 (the upper levels of the inter- 
combination components) begin to dominate over the ra- 
diative transition from ^5*1 to the ground level. Thus, on 
the specified grid of parameters no <C ^e,cnY, i.e., the in- 
tensity redistribution of the forbidden and intercombination 
lines for iron may be neglected. Note that this conclusion 
is also valid for the nickel triplet, because rie^crit = ^ ocZ 



( Mewe. Schrijveij|T978a ). where A/^ is the rate of radiative 
transitions to the ground level and C// is the rate of colli- 
sional excitations to the upper levels of the intercombina- 
tion lines. 

The results of our analysis of the simulated spectra in 
the Chandra HEG resolution are presented in Table |3] 
Since FWHMch dXE^'^l ke V exceeds the expect ed line 



FWHMo = V^7/3esinx E, ([Marshall et al.l l2002h (x 
65° is the angle between the jet axis and the observer's 
direction) for 6 0.01 rad (see Fig. O, we fitted the 
spectrum convolved with the model response function by 
three Gaussians of the same FWHM but with decoupled 
centroids {Ef^Ei^Ey in the rest frame of the jet) and ampli- 
tudes/ /, r (without scattering {f/r)^ = 0.28, (//r)^ = 0.27 
(see TableO. In what follows, the tables of results give the 
line width W minus the instrumental broadening (i.e., W = 

^JfWHM^ - FWHMl^ and W = ^JEWHM^ - EWHM\^ 
in the cases of analysis in the Chandra HEG and Astro-H 
resolutions, respectively). 

The attenuation coefficient (1 - C) was calculated for the 
triplet as a single line with a centroid Ec (the first mo- 
ment of the photon number distribution in energy), i.e., its 
boundaries were determined with respect to the common 
broad wings (see Table [3]). 

Since the spectra convolved with the Astro-H model re- 
sponse function cannot be described in such a way, we 
provide only some of the derived line profiles (Fig. O, 
which illustrate the characteristic features of a particular 
scattering regime. The energy E in the observer's frame 
of reference is along the horizontal axis; the fraction of all 
line photons with energies in the range from E - ISE /2 to 
E + /S.E /2, where /S.E <C FWHM is the spectral bin size, is 
along the vertical axis. 

As has been noted above, the triplet components at 9 > 



0.01 rad actively interact via the optical depths, what dis- 
torts noticeably the line shape. For 6 = 0.02 rad the trans- 
verse optical depth of the jet for Thomson scattering tt 
turns out to be less than a critical value r^r ~ 0.19 (see the 
"Radiative Transfer in Lines" Section) for the resonance 
line photons at ro > 4 x 10^^ cm; therefore, the scatterings 
by electrons cause all triplet lines to be suppressed to ap- 
proximately the same degree. The fraction of the photons 
that left the line C reaches 30% at ro = 4 x lO^^cm. Note that 
for 6 = 0.02 rad the quantities Ef^Et^Er and/ /, r found in 
our analysis lose their physical linkage to the forbidden, in- 
tercombination, and resonance lines, respectively, and ba- 
sically have only a descriptive significance. 

For 6 = 0.01 rad the interaction between the components 
still remains significant, but the "extension" of the optical 
depth to the low-energy region for the resonance line pho- 
tons becomes its main manifestation. In this case, tt turns 
out to be greater than r^r ^ 0.15 for the resonance pho- 
tons already at ro = 8 x 10^^ cm. This enhances the escape 
into the wings ( and increases the relative intensities of the 
forbidden (f/r) and intercombination (i/r) lines (see Table 
O. Thus, the fraction of the photons that left the line turns 
out to be significant, but, nevertheless, ( < 40% for the 
entire triplet due to the photons emitted inside the surface 
layer (see the "Radiative Transfer in Lines" Section) and 
the contribution from the forbidden and intercombination 
components. Another effect is the increase in line width W 
as a result of multiple scatterings. 

For 6 = 0.003 rad the interaction between the compo- 
nents no longer plays a significant role, while tt turns out 
to be greater than r^r ^ 0.1 1 for the resonance photons al- 
ready at neo> I X 10^"^ cm~^. Therefore, the resonance line 
is suppressed noticeably in comparison with the forbidden 
and intercombination lines. The existence of a transverse 
velocity gradient has an effect on the strong asymmetry of 
the resonance line profile, in agreement with the predic- 
tions made in the "Radiative Transfer in Lines" Section. 
The additional broadening also becomes more pronounced. 

In the cylindrical case of 9 = 0.0001 rad, the transverse 
velocity gradient is negligible and almost all of the reso- 
nance line photons produced deep in the jet are ultimately 
scattered by an electron and fall into the broad line wings. 
Thus, in a sense, a pure case is realized, i.e., only the pho- 
tons emitted in the jet surface layer remain at the resonance 
line center. As a result, the entire triplet noticeably loses 
in intensity, while the relative contribution of the forbid- 
den and intercombination lines increases. In this case, the 
width of the lines exceeds that of the unscattered lines at 
6 = 0.003 rad and is less than that of the scattered lines at 
6 = 0.01 rad by only a factor of 1.5. Given the additional 
possibilities for broadening (see Section 5.3), this indicates 
that 6 is difficult to measure accurately based on Chandra 
HETGS observations of the helium-like iron triplet. In 
turn, Astro-H is ideally suited for studying the fine struc- 
ture of this line and diagnosing the jet parameters from it. 
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Table 2. The grid of input parameters and the corresponding values of o; at Tq = 20 keV as well as the luminosities (L„ = L/10'^) and equivalent 
widths (without scattering) of the lines from Table 1 for which detailed simulation results are presented (see the "Results" Section). The use of the 
quasi-cylindrical calibration (|9) when passing from tq = 8 x 10^*^ cm to tq = 4 x lO^^cm for 6 = 0.01 and 6 = 0.02 rad is marked (*) , because a ~ 1 
in this case 









Fe XXV Kc, 


Fe XXVI Lyc, 


Ni XXVII Kc, 


Fe XXV Kb 


SXV Kc, 


Si XIII Kc, 




10^4 


a 


L33, 


EW 


L33 


EW 


L32 


EW 


L32 


EW 


L33 


EW 


L33 


EW 


cm 


cm~^ 




erg/s 


eV 


erg/s 


eV 


erg/s 


eV 


erg/s 


eV 


erg/s 


eV 


erg/s 


eV 


Quasi-adiabatic models 














6 


= 0.02 rad 
















10.0 


0.16 


0.10 


8.7 


227 


5.00 


164 


5.59 


21 


9.11 


34 


0.99 


13.3 


1.48 


15.4 


6.4 


0.31 


0.13 


8.4 


230 


4.73 


163 


5.34 


21 


8.74 


35 


0.91 


13.7 


1.44 


15.8 


3.6 


0.75 


0.17 


8.4 


234 


4.63 


163 


5.34 


22 


8.7 


35 


0.85 


14.2 


1.44 


15.9 


1.6 


2.5 


0.25 


7.4 


239 


4.00 


162 


4.66 


22 


1.1 


36 


0.65 


14.3 


1.23 


15.4 


0.8 


lO.Ol 


0.51 


11.1 


242 


5.88 


162 


6.91 


22 


11.47 


37 


0.98 


13.6 


1.61 


13.7 














6 


= 0.01 rad 
















5.0 


0.32 


0.20 


7.8 


238 


4.29 


163 


4.96 


22 


8.17 


36 


1.04 


14.3 


1.34 


15.7 


3.2 


0.62 


0.25 


7.3 


239 


3.94 


162 


4.58 


22 


7.57 


36 


1.03 


14.3 


1.21 


15.4 


1.8 


1.5 


0.34 


6.9 


241 


3.68 


162 


4.3 


22 


7.13 


36 


1.05 


14.1 


1.09 


14.8 


0.8 


5.0 


0.51 


5.6 


242 


2.95 


162 


3.46 


22 


5.73 


37 


0.93 


13.6 


0.80 


13.7 


0.4 


20.01 


1.02 


7.0 


240 


3.76 


162 


4.38 


22 


7.24 


36 


1.30 


12.4 


0.89 


12.1 


Quasi-cylindrical models 
















= 0.003 rad 
















2.0 


1.0 


8.69 


9.0 


241 


4.84 


162 


5.62 


22 


9.3 


36 


0.99 


12.6 


1.17 


12.4 


1.0 


4.0 


17.38 


10.1 


237 


5.54 
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6.4 


22 


10.34 


36 


1.04 


11.9 


1.22 


11.3 


0.5 


16.0 


34.77 


10.8 
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6.03 
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6.89 


22 


11.31 


35 


1.08 


11.4 


1.24 


10.8 














e = 


0.0001 rad 
















2.0 


1.0 


26.07 


11.5 


231 


6.59 
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7.46 


21 


12.22 


35 


1.12 


10.9 


1.28 


10.3 


1.0 


4.0 


52.15 


11.5 


231 


6.52 
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7.35 


21 


12.03 


35 


1.10 


10.8 


1.26 


10.3 


0.5 


16.0 


104.30 


11.5 


231 


6.52 
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7.37 


21 


12.06 


35 


1.10 


10.8 


1.26 


10.4 
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Table 3. Analysis of the simulation results for the FeXXV Ka triplet in the Chandra resolution for the approaching jet with Tq = 20 keV at a phase 
corresponding to Zb = -0.078. The widths of the unscattered lines are Wq = 37.8 eV at 6 = 0.02 rad, Wq = 24.0 eV at 8 = 0.01 rad, Wq = 13.3 eV at 
e = 0.003 rad, and = 11-9 eV at 8 = 0.0001 rad 







1-C 


Ec, keV 


W,eV 


f/r 


i/r 


E/, keV 


E„ keV 


E„ keV 


^0 


10^^ cm 








6= 


=0.02 rad 










5.0 


0.95 


6.6835 


38.2 
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0.630 


6.6336 


6.6736 


6.7048 
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0.4 


0.69 
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6.6721 
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1.8 
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6.6831 


25.9 


0.298 


0.274 


6.6378 


6.6721 


6.6989 




0.8 


0.81 


6.6829 


26.1 


0.307 


0.275 


6.6384 


6.6721 


6.6989 




0.4 


0.62 


6.6826 


26.2 


0.338 


0.281 


6.6397 


6.6730 


6.6995 












e= 


=0.003 rad 










1 


0.86 


6.6824 


15.9 
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0.300 


6.6372 


6.6731 


6.6985 




4 


0.71 


6.6813 


16.8 
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0.333 


6.6376 


6.6741 


6.6984 




16 


0.51 


6.6800 


17.3 
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0.408 


6.6381 


6.6758 


6.6990 




10^4 cm-^ 








e= 


0.0001 rad 










1 


0.80 


6.6820 


14.7 


0.338 


0.334 


6.6370 


6.6733 


6.6994 




4 


0.64 


6.6803 


15.4 


0.406 


0.392 


6.6372 
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6.6996 




16 


0.45 


6.6781 


15.6 


0.515 


0.486 


6.6373 


6.6758 


6.6999 




Fig. 5. Results of our simulations for the FeXXV triplet for the approaching jet with Tq = 20 keV at a phase corresponding to = -0.078. (a) 
8 = 0.01; the upper panel shows the spectrum convolved with the Astro-H model response function; the lower panel show the spectrum convolved 
with the Chandra HEG model response function; the dotted curve for rro = 4 x 10^^ cm, the curve with short dashes for tq = 8 x 10^^ cm, the curve 
with long dashes for tq = 1.8 x 10^^ cm, the dash-dotted curve for ro = 3.2 x 10^^ cm, and the solid curve indicates the unscattered line profile (the 
same for all tq). (b) 8 = 0.003 rad; the upper panel shows the spectrum convolved with the Astro-H model response function; the lower panel shows 
the spectrum convolved with the Chandra HEG model response function; the dotted curve for Ugo = 1.6 x 10^^ cm~^, the curve with short dashes for 
«eO = 4 X lO^^cm"^, the dash-dotted curve for Ugo = 1.0 x 10^^ cm~^, and the solid curve indicates the unscattered line profile (the same for all Ueo) 
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5.1.2 FeXXVILya 

Since the FeXXVI Lyo, doublet corresponds to the tran- 
sitions in hydrogenlike iron, both emissivity and optical 
depth of the jet plasma in these lines have a peak at slightly 
higher temperatures than those in the helium-like iron 
triplet lines (see Fig. |4]). This allows the Fe XXVI LyJ Fe 
XXV Kq, ratio to be used to dia gnose the t emperature of 
the hottest parts of the jet (IKotani et alJI 19961) . 

In this case, the oscillator strength / for the Lyo, transi- 
tions is considerably (approximately by a factor of 3 for the 
ls-2p transition and approximately by a factor of 6 

for the \s-2p CP1/2) transition; see Table [TJ smaller than 
that for the resonance component of the triplet, whence, 
at first glance, the doublet lines can be assumed to be less 
subjected to the scattering effects inside the jet. However, 
the Doppler energy shift due to the transverse jet plasma ve- 
locity component at 6 ~ 0.01 rad allows the photons of the 
higher-energy doublet component to be efficiently scattered 
on the transition corresponding to the low-energy compo- 
nent. As a result, the scattered line for 6 = 0.01 rad has 
a completely different profile compared to the initial one 
(see Fig. 6); for its description in the Astro-H resolution, 
we used a model consisting of two Gaussians of the same 
FWHM with amplitudes Al and A2 and centroids shifted 
by 5zi and Szi relative to the initial positions for the low- 
and high-energy components, respectively (Table O. 

When analyzing the spectra in the Chandra HEG reso- 
lution, we used a model of a single Gaussian with FWHM 
with a centroid corresponding to the weighted mean energy 
of the unscattered doublet (Tabled. The attenuation coeffi- 
cient (1 - is given for the doublet as a single line. Typical 
scattered line profiles are shown in Fig. [6] 

Although for 6 = 0.02 rad and 6 = 0.01 rad tt is less 
than Tcr for the photons of the high-energy component at 
all ro, except for ro = 4 x 10^^ cm, the total doublet inten- 
sity decreases considerably. This is probably a result of the 
above-mentioned extension of the optical depth to the red 
region for the photons of the high-energy component. This 
effect is clearly illustrated by the results of our analysis of 
the scattered doublet profiles in the Astro-H resolution (Ta- 
ble O. The increase of A1/A2 with decreasing ro allows the 
"flattening" of the doublet profile to be judged, while the 
centroid shifts 5zi and Szi point to the absence of a phys- 
ical linkage of the fitting Gaussians to the doublet compo- 
nents. In addition, the effective line broadening turns out 
to be very significant (Table |4]). 

For 6 = 0.003 rad and 9 = 0.0001 rad, in the absence 
of any interaction between the components, the scattering 
pattern is completely identical to that described in detail for 
the Ko, triplet. In this case, the scattered line widthW from 
our analysis in the Chandra HEG resolution turns out to be 
approximately equal to the unscattered line width Wo for 
e = 0.01 rad. 

Thus, owing to the interaction between the components. 



the influence of scattering effects turns out to be significant 
(primarily on the line profiles) even at large 6 almost on the 
entire grid of parameters. This makes the fine structure of 
the Ljo, doublet profile a very sensitive tool for diagnosing 
the hottest parts of the jet, in particular, based on Astro-H 
observations. 

Table 5. Analysis of the simulation results for the FeXXVI Ly^ doublet 
in the Astro-H resolution for the approaching jet with Tq = 20 keV corre- 
sponding to = -0.078. 



e = 0.01 rad 



ro,10^^ cm 


A1/A2 


W, eV 


fei,io-4 


&2,10-^ 


5.0 


0.66 


26.18 


1.46 


-1.69 


3.2 


0.69 


25.73 


1.34 


-2.17 


1.8 


0.73 


25.57 


1.59 


-2.66 


0.8 


0.78 


24.96 


2.20 


-3.52 


0.4 


0.80 


24.78 


1.48 


-5.20 



5.1.3 Ni XXVII Ka +Fe XXV 

The currently available observations near the nickel Kq, 
triplet point to an excess of radiation compared to that ex- 
pected from the FeXXV (6.7 keV) line intensity within 
the standard model with solar elemental abundances (Ta- 
ble. [2]). This is interpreted i n terms of a nick e l over - 
abundance i n the jet plasma (Brinkmann et al. ('2005"), 
Medvedev, Fabrikal (|2010)). Since present-day instruments 
do not aUow NiXXVII K^, and FeXXV to be reliably 
resolved, we performed joint simulations of these lines, al- 
though a noticeable overlap takes place only for 6 = 0.02 
rad (see Fig. (T]). 

We do not provide the results of our analysis of the simu- 
lated spectra in the Chandra HEG resolution due to the low 
sensitivity of this instrument near E^'^^.S keV. 

When analyzing the spectra in the Astro-H resolution, 
we used a model consisting of five Gaussians of the same 
FWHM with decoupled amplitudes but fixed centroids. 
The attenuation coefficients {l-Q are determined sepa- 
rately for the nickel triplet and K/3 of helium-like iron (see 
TableO. 

Whereas tt is much less than Tcr for the photons of the 
nickel resonance line on the entire grid of quasiadiabatic 
model parameters, Tcr is not so large for the photons of the 
iron line and even turns out to be less than fr at 6 = 0.01 
rad and ro = 4 x 10^^ cm. In particular, this is reflected 
in different dependences of C on ro (see Table refnifeah). 
As a result, the F(Ni)/F(Fe) ratio slightly changes, but this 
change takes place only at ro = 4 x 10^^ cm and does not 
exceed 10% in magnitude. Thus, the main effect is a gen- 
eral decrease in intensity due to Thomson scattering almost 
on the entire grid of quasi-adiabatic model parameters. 

For the quasi-cylindrical models, the overall pattern of 
radiative transfer is very similar to the case of the FeXXV 
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Table 4. Analysis of the simulation results for the FeXXVI Lja doublet in the Chandra resolution for the approaching jet with Tq = 20 keV correspond- 
ing to Zb = -0.078. The widths of the unscattered lines are Wq = 54 J eV at 8 = 0.02 rad , Wq = 35.9 eV at 8 = 0.01 rad , Wq = 29.6 eV at 8 = 0.003 rad 
, Wo = 29.5 eV at 8 = 0.0001 rad . 



1 r\ 1 1 

ro, 10 cm 
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6=0.02 rad 
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=0.01 rad 


W, eV 
9=0.02 rad 
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:0.01 rad 


5.0 


0.93 




0.93 


59.2 




40.3 
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0.91 
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61.1 




40.9 


1.8 


0.88 




0.86 


62.3 




41.8 


0.8 


0.82 




0.78 


64.5 




43.7 


0.4 


0.65 




0.58 


68.3 




45.3 


Heo^ 10^4 cm"-^ 


6 = 3 X 10-^ rad 
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= lO-^rad 


e = 3 X 10-^ rad 


e = 


10-4 rad 


1.0 


0.86 




0.82 


30.9 




31.1 


4.0 


0.71 




0.66 


31.8 




32.0 


16.0 


0.49 




0.45 


32.9 




33.0 



Ka triplet, given that for the resonance line of the nickel 
triplet and Kf^ of iron tt > Tcr only for a > 30 and a > 10, 
respectively (see Table [2]) . 

Thus, the intensity "transfer" between the iron Kjs line 
and the nickel Ko, triplet through resonant scattering turns 
out to be fairly likely. At the same time, the intensity of 
the entire set of lines does not undergo such a dramatic 
decrease as the iron Koc triplet. However, this cannot com- 
pletely explain the observed excess of radiation in this re- 
gion relative to that expected from the FeXXV Ko, (6.7 
keV) line under the assumption of solar elemental abun- 
dances. 

Astro-H observations allow the nickel triplet lines and 
Kj3 of iron to be reliably resolved, thereby making it possi- 
ble to determine the nickel abundance relative to iron in the 
jets of SS 433 in situ. 

5.1.4 SXVKa 

The overall pattern typical of the Kq, triplets of helium- 
like ions at low energies (we will be interested primarily in 
sulfur and silicon; see Table [2]) differs from that for high- 
energy triplets (e.g., of iron) in several respects. First, the 
triplet components have a large relative spectral separation 
(AE / AEd)', as a result, the overlap and, hence, the inter- 
action take place only for large 6 ~ 0.02 rad (see Fig. [8^). 
Nevertheless, we provide the attenuation coefficients l-( 
for the entire triplet (see Table IT). Second, Ue^crit oc Z^^ 
turns out to be 10^^ cm (an accurate calculation for the 
corresponding temperature range was performed using the 
Cloudy 08.00 code). Therefore, for the fine structure of the 
triplets to be simulated in detail, the intensity redistribution 
of the forbidden and intercombination lines as a result of 
the collisional excitations from ^5"! (the upper level of the 
forbidden component) to ^Po,i,2 (the upper levels of the in- 
tercombination comp onents) should be taken into account 
( Porquet et al.l |201Q|) . This redistribution is described in 



terms of their intensity ratio R(ne) = z/(x-l-y) = 



l+ne/ne,crit(T)' 



where ro is the value of this ratio in the low-density limit 
(from APEC), and the dependence Ue^crit (T) is obtained by 
fitting the Cloudy computation results. 

Since Astro-H is inferior to Chandra HETGS in its spec- 
tral characteristics near the brightest lines of hydrogen- and 
helium-like sulfur and silicon, we omit the results of our 
analysis of the simulated spectra in the Astro-H resolu- 
tion. For our analysis in the Chandra HEG resolution (Fig. 
[8^) we used, where possible, a model consisting of three 
Gaussians of the same FWHMwith fixed centroids but with 
decoupled amplitudes w, x + y, and z for the resonance, 
intercombination, and forbidden components, respectively. 
Since the triplet components have different characteristics 
with respect to the scattering effects, one might expect 
these effects to affect noticeably the ratio R = z/(x + y) 
and the ratio G(T) = (z + x + y)/w of the total intensity 
of the forbidden and intercombination lines to the intensity 
of the resonance line, which is used t o diag nose the plasma 
temperature (see, e.g. JPorquet et al.l ( 2010l) ). Given the im- 
portance of these ratios from the viewpoint of the physical 
interpretation of observational data, the summary table of 
results provides precisely the values of R and G (see Table 

m. 

For 6 = 0.02 rad the interaction between the triplet 
components resembles the situation with the Ka triplet of 
helium-like iron at O = 0.01 rad, but complicated by the 
intensity redistribution of the forbidden and intercombina- 
tion components. Therefore, we do not present the results 
of a detailed analysis for 6 = 0.02 rad, bearing in mind the 
impossibility of their clear and unambiguous interpretation. 

For 6 = 0.01 rad, as a result of the rapid decrease in 
He oc with increasing ^, the effective optical depth for 
scattering by free electrons tt turns out to be greater than 
Tcr for the resonance photons only at ro = 4 x 10^^ cm, de- 
spite the large optical depth of the jet for resonant scattering 
in the allowed S XV triplet component (see Fig. |4]). There- 
fore, the overall pattern is similar to that considered for the 
Ni XXVII Ka triplet, i.e., the influence of scattering effects 
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Fig. 7. Results of our simulations for the Ni XXVII and Fe XXV 
Kj^ for the approaching jet with Tq = 20 at a phase corresponding to Zb = 
-0.078. (a) The spectrum convolved with the Astro-H model response 
function for O = 0.02 rad; (b) the spectrum convolved with the Astro- 
H model response function for O = 0.01 rad; the dotted curve for tq = 
4 X 10^^ cm, the curve with short dashes for tq = 8 x 10^^ cm, the dash- 
dotted curve for tq = 1.8 x 10^^ cm, and the solid curve indicates the 
unscattered line profile (the same for all tq); (c) the spectrum convolved 
with the Astro-H model response function for 6 = 0.003 rad; the dotted 
curve for = 1-6 x 10^^ cm the curve with short dashes for Ueo = 
4 X 10^^ cm the dash- doted curve for Ueo = 1.0 x 10^^ cm and the 
solid curve indicates the unscattered line profile (the same for all Ueo). 



is small, the ratios R and G differ from their unscattered 
values Ro and Go only slightly (see Table (T]). 

In the case of quasi-cylindrical models, the overall pat- 
tern of the spectra corresponds to that for the Fe XXV 
triplet (here, rie changes little). As a result, the ratios R 
and G as well as the effective line width change noticeably. 
However, even in the case of an increase by a factor of 1.5- 
2, the line width for O ^ 0.003 rad remains considerably 
smaller than that at 6 = 0.01 (as distinct from the analo- 
gous situation for the FeXXV triplet). 

5.7.5 Si XIII Ka 

The general scheme for analyzing and presenting the re- 
sults for the Ka triplet of helium-like silicon is completely 
identical to the scheme described in detail for the sulfur 
triplet in the preceding section. Therefore, here we omit 
the determination of the quantities given in Table O Simi- 
larly, the interpretation of the results of our analysis closely 
coincides with that for the sulfur triplet. 

Nevertheless, an important distinctive feature of the sil- 
icon triplet is the possible influence of satellites on mea- 
sured characteristics (e.g., R and G). We will defer a de- 



tailed discussion of this question until Section [7. 3. II while 
for now note that the diagnostics of plasma parameters 
based on the ratios R and G can be strongly affected by 
the scattering effects. For example, an increase in these ra- 
tios (see Tables [7] and (8]) as a result of scattering relative to 
the v alues correspo nding to the optically thin approxima- 
tion ('Porque t et al.ll201 0) must lead to an underestimation 
of the density and temperature from observations. 

Given the excellent spectral resolution of Chandra 
HETGS in this energy range, a detailed analysis of the line 
profiles corrected for the effects considered can provide in- 
formation about the low-temperature part of the jet. 

5.2 Broad Line Wings 

Broad line wings are formed at a large number of pho- 
tons that left the line due to their scattering by a free elec- 
tron. The shape of these wings is well described by the 
single Compton scattering kernel with the application of 
relativistic corrections in the ca se of an isotrop ic radiation 
field defined by Eq. (19) from ISazonov. Sunyaev. (20001 
In this case, the only normalization is the quantity (, i.e., 
the fraction of the photons that left the line, while the only 
parameter in fitting is the effective temperature of the scat- 
tering electrons Tg. This temperature reflects the contri- 
bution from various parts of the jet to the scattered radia- 
tion and, hence, allows the temperature and density profiles 
along the jet to be judged. For example, on average, higher- 
temperature fits correspond to the class of quasi-adiabatic 
models than to the class of quasi-cylindrical ones (see Fig. 
[9]). This is due to a rapid decrease in density against the 
background of a slow change in temperature for the former 
and the directly opposite situation for the latter. 

This single Compton scattering kernel has the following 
characteristic features: first, the presence of a cusp near 
E = Eq and, second, a noticeable asymmetry- the "right" 
(high-energy) wing is much broader than the "left" (low- 
energy) one. In addition, there is a probability of photon 
rescattering by an electron. However, at fr ~ 0.1 the influ- 
ence of rescatterings is negligible, with the possible excep- 
tion of the far edges of the wings (see Fig. O. 

Since, in general, the broad wings of scattered lines are 
always present in the spectrum(the question is only in what 
quantities), an accurate measurement of their shape can 
be used as a universal (though model-dependent) tool for 
determining the jet parameters that is already accessible 
for present-day instruments. Undoubtedly, the contribution 
from the wings of neighboring lines slightly complicates 
the problem, but, at the same time, independent measure- 
ments in the high- and low-energy parts of the spectrum 
can reveal particular local features in the density and tem- 
perature distributions. 

At the same time, it should be remembered that the con- 
tinuum itself is also subject to scattering inside the jet, 
which is especially important for the photons in resonance 
with the electron transitions in ions. Nevertheless, our sim- 
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ulations with the addition of an appropriate wide contin- 
uum component showed that the continuum "subsidence" 
under the resonance hnes for Hnes with significant equiva- 
lent widths affects only slightly the observed parameters of 
these lines (^5% for the intensity). Since a detailed dis- 
cussion of the scattering effects on the continuum is beyond 
the scope of this paper, we will restrict ourselves here only 
to the remark made and will consider a number of effects 
affecting directly the lines. 



0=10--' rad To = 20 keV 
n =1x1014 cm-3 




9 10 



Energy, keV 




Table 6. Analysis of the simulation results for the Ni XXVII +Fe 
XXV in the Astro-H resolution for the approaching jet with Tq = 20 
keV corresponding to = -0.078. The widths of the unscattered lines are 
Wo = 30.3 eV at e = 0.01 rad , Wq = 10.8 eV at 8 = 0.003 rad , Wq = 7.65 
eV at e = 0.0001 rad . 





l-C(Ni) 


l-C(Fe) F(Ni)/F(Fe) 


W, eV 


ro,10^^ cm 




6=0.01 rad 




5.0 


0.94 


0.94 0.61 


30.3 


3.2 


0.92 


0.92 0.61 


30.6 


1.8 


0.89 


0.89 0.61 


30.9 


0.8 


0.84 


0.83 0.61 


31.3 


0.4 


0.71 


0.65 0.65 


31.3 


Heo^ 10^^ cm ~^ 




6=0.003 rad 




1 


0.90 


0.88 0.62 


13.2 


4 


0.80 


0.74 0.65 


14.7 


16 


0.63 


0.53 0.72 


16.2 


Heo^ 10^^ cm ~^ 




6=0.0001 rad 




1 


0.87 


0.83 0.64 


11.2 


4 


0.77 


0.68 0.70 


12.6 


16 


0.60 


0.48 0.77 


14.1 



Energy, keV 



Fig. 9. Fitting the broad wings of the scattered FeXXV Ka line in the 
jet's frame of reference by the single Compton electron scattering ker- 
nel from Sazonov, Sunyaev ( 2000) for one of the quasi-cylindrical models 
(top) and one of the quasi-adiabatic models (bottom). The photon scatter- 
ing probability for a line with the central energy Eq = 6.1 keV and a cloud 
of electrons with temperature Te is along the vertical axis. 
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Fig. 6. Results of our simulations for the FeXXVI Lya doublet for the approaching jet with Tq = 20 keV at a phase corresponding to zt = -0.078. 
(a) 6 = 0.01; the upper panel shows the spectrum convolved with the Astro-H model response function; the lower panel show the spectrum convolved 
with the Chandra HEG model response function; the dotted curve for rro = 4 x 10^^ cm, the curve with short dashes for tq = 8 x 10^^ cm, the curve 
with long dashes for tq = 1.8 x 10^^ cm, the dash-dotted curve for tq = 3.2 x 10^^ cm, and the solid curve indicates the unscattered line profile (the 
same for all ro). (b) O = 0.003 rad; the upper panel shows the spectrum convolved with the Astro-H model response function; the lower panel shows 
the spectrum convolved with the Chandra HEG model response function; the dotted curve for Ueo = 1.6 x 10^^ cm~^, the curve with short dashes for 
neo = 4x lO^'^cm"^, the dash-dotted curve for Ueo = 1.0 x 10^^ cm~^, and the solid curve indicates the unscattered line profile (the same for all Ueo) 



Table 7. Analysis of the simulation results for the SXV triplet in the Chandra resolution for the approaching jet with Tq = 20 keV at a phase 
corresponding to = -0.078 







1-C 


R(Ro) 


G(Go) 


W (Wo) , eV 


^0 


,10' ' cm 




e 


= 0.01 rad 






5.0 


0.97 


1.98 ( 1.92) 


0.65 ( 0.64 ) 


8.92 ( 8.77 ) 




3.2 


0.96 


1.98 ( 1.92) 


0.65 ( 0.64 ) 


8.96 ( 8.77 ) 




1.8 


0.93 


1.91 ( 1.92) 


0.65 ( 0.64 ) 


9.06 ( 8.77 ) 




0.8 


0.88 


1.58 ( 1.62) 


0.65 ( 0.64 ) 


9.09 ( 8.77 ) 




0.4 


0.71 


0.71 (0.73) 


0.67 ( 0.64 ) 


8.97 ( 8.77 ) 




10^4 cm-^ 




e = 


= 0.003 rad 






1 


0.90 


1.87 ( 1.87) 


0.64 ( 0.59 ) 


3.57 (2.85) 




4 


0.75 


1.39 ( 1.34) 


0.69 ( 0.59 ) 


3.83 (2.85) 




16 


0.55 


0.66 ( 0.57 ) 


0.76 ( 0.59 ) 


4.16(2.85) 




10^^ cm-^ 




e = 


0.0001 rad 






1 


0.80 


1.75 ( 1.71) 


0.71 (0.59) 


2.93 ( 1.73) 




4 


0.65 


1.26 ( 1.14) 


0.81 (0.59) 


3.25 ( 1.73) 




16 


0.46 


0.62 ( 0.49 ) 


0.90 ( 0.59 ) 


3.52 ( 1.73) 
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Fig. 8. Results of our simulations for the SXV Kqc triplet (a) and the SiXIII Kql triplet (b) for the approaching jet with Tq = 20 keV at a phase 
corresponding to Zb = -0.078. The upper panel: B = 0.02 rad; the spectrum convolved with the Chandra HEG model response function; the dotted curve 
for ro = 4 X 10^^ cm, the curve with short dashes for tq = 8 x lO^^cm, the dash-dotted curve for ro = 1.8 x 10^^ cm, 1, 2, 3 represent the unscattered 
line profile for tq = 4 x 10^^ cm, tq = 8 x 10^^ cm, and tq = 1.8 x 10^^ cm, respectively. The middle panel: B = 0.01 rad; the spectrum convolved with 
the Chandra HEG model response function; the dotted curve for tq = 4 x 10^^ cm, the curve with short dashes for tq = 8 x 10^^ cm, the dash-dotted 
curve for tq = 1.8 x 10^^ cm, 1, 2, 3 represent the unscattered line profile for tq = 4 x 10^^ cm, tq = 8 x 10^^ cm, and ro = 1.8 x 10^^ cm, respectively. 
The lower panel: B = 0.003 rad; the spectrum convolved with the Chandra HEG model response function; the dotted curve for neo = 1 x 10^^ cm 
the curve with short dashes for neo=4x 10^^ cm the dash-dotted curve for Ugo = 1.6 x 10^^ cm 1, 2, 3 represent the unscattered line profile for 
Ueo = 1.6 X 10^^ cm ~^,neo=4x 10^^ cm and neo = lx 10^^ cm respectively. 



Table 8. Analysis of the simulation results for the SiXIII Kct triplet in the Chandra resolution for the approaching jet with Tq = 20 keV at a phase 
corresponding to Zb = -0.078 





1-C 


R(Ro) 


G(Go) 


W (Wo) , eV 


ro,10^^ cm 




e 


= 0.01 rad 




5.0 


0.98 


2.33 (2.33) 


0.720(0.713) 


6.90 ( 6.76 ) 


3.2 


0.96 


2.25 ( 2.27 ) 


0.716(0.709) 


6.92 ( 6.75 ) 


1.8 


0.94 


1.95 ( 1.98) 


0.709 ( 0.699 ) 


6.91 (6.71 ) 


0.8 


0.89 


1.09 ( 1.10) 


0.706(0.691) 


6.93 ( 6.67 ) 


0.4 


0.72 


0.26 ( 0.26 ) 


0.708(0.671 ) 


6.68 ( 6.56 ) 


Heo, 10^^ cm~^ 




e 


= 0.003 rad 




1 


0.90 


1.79(1.76) 


0.70 (0.64) 


2.63 (2.08) 


4 


0.76 


0.73 (0.67) 


0.74 (0.64) 


2.80 (2.08) 


16 


0.55 


0.22 (0.17) 


0.79 (0.64) 


2.96 (2.08) 


Heo, 10^^ cm~^ 




e = 


= 0.0001 rad 




1 


0.78 


1.38(1.24) 


0.81 (0.63) 


2.26(1.15) 


4 


0.63 


0.61 (0.48) 


0.89 (0.63) 


2.51 (1.15) 


16 


0.45 


0.21 (0.14) 


0.91 (0.63) 


2.68(1.15) 
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5. 3 The Influence of Other Effects 

Apart from the scattering inside the jet, there are also 
several effects capable of affecting the observational line 
characteristics (intensity, width, etc.) and the line ratios 
used to diagnose the plasma parameters. In this Section, 
we estimate the possible effects of satellites of bright lines, 
nutations, and photoionization and photoexcitation. 

5.3.1 Satellites and the Overlap of Close Lines 

On the one hand, the presence of satellites near bright 
lines slightly complicates the analysis and interpretation of 
the spectra, but, on the other hand, investigating the satellite 
parameters (provided an appropriate spectral resolution of 
the instrument) is a powerful tool for diagnosin g the param- 



eters o f the emitting medium per se (see, e.g. JPorquet et al 
(I20TqI) ). 



For the triplets of helium-like ions, the bulk of the satel- 
lite intensity integrated over the entire jet is concentrated in 
about 25 lines inside the triplet, i.e., between the forbidden 
and resonance lines. In this case, the satellite intensity rela- 
tive to the resonance line slightly increases with increasing 
charge of the ion nucleus Z. For example, the contribution 
of the satellites in the jet model under consideration for sul- 
fur (Z = 16) is about 12% of the intensity of the Ko, triplet 
of the helium-like ion and about 20% for the same triplet 
of iron (Z = 26). At the same time, there can be cases 
where the energies of the bright lines corresponding to the 
transitions in various ions of the same or different elements 
are close. In particular, such a situation takes place for the 
forbidden line of the helium-like silicon Kq^ triplet coinci- 
dent with the MgXII Ly^ doublet. As a result, not only the 
normalization of the ratios R and G (see Section IF. 1.41 ) but 
also the shape of the temperature dependence G(T) itself 
changes (see Fig. [TOh). 

In the case of intrinsically broad lines, the satellites can- 
not be resolved in principle, but their contribution to the 
total intensity of the blend and its width can be substantial. 
For example, for the triplet of helium-like iron, the effec- 
tive increase in the width W of the blend components (for 
a definition of W, see Section 15.1.11) due to the contribu- 
tion of satellites turns out to be /SW > 1 eV, i.e., > 5% at 
6 = 0.01rad (see Fig.[TOb). 

5.3.2 Nutations 

The nutation and, possibly, jitter of the jets 



(iKubota et al.l 1201 Ob are other factors contributing to the 
broadening of the observed lines. To estimate the nuta- 
tion b roadening, we used the mode l of ephem e rides ( Katzl 



1987h with the parameters found bv lOies et all d 2002k that 



defines z for each of the jets at an arbitrary instant of time. 
The influence on the line width can be described in terms of 
the total change Az in some time interval (^o,^ +^oZ?^) short 
compared to the nutation period (6.28 days), where is the 



Table 9. Nutation broadening of the Fe XXV (6.70 keV) and Si XIII (1.86 
keV) lines for phases correspondi ng to the Chan dra observations, a repor t 
on which is contained in Marshal l et al.l ( |2002|) [1], iNamiki et ( l2003h 
[2], and L opez et all (t2006>) [3]. 



Report 


Zb 


^obs-) ks 


Az,10-4 


5Eq, eV 
Fe Si 


[1] 


-0.078 


32.0 


5 


8 2.2 


[2] 


0.046 


21.3 


9 


11 3.1 


[3] 


0.014 


25.7 


10 


13 3.6 



starting time of observations and tobs is their total duration. 
For example, for tobs = 20 ks, the maximum value of Az is 
about 0.0024, which is in good agreement with the mea- 
sure d maximu m nutation variability dz/dt ^ 0.004 in 0.33 
day i Fabrikat2004 ). In this case, as long as the total time of 
observations is much shorter than the nutation period (6.28 
days), Az scales linearly in tobs with a good accuracy. The 
change in the energy Eq of the line centroid due to nutation 
can be estimated as SEq = Az. The full width at half 



maximum of the line with nutation is then 



FWHM„ 



FWHM^ + ^^SE^, 



(24) 



where FWHM is the full width at half maximum without 
nutations. The values of Az and SEq for the Fe XXV (6.7 
keV) and Si XIII (1.86 keV) lines of the "blue" jet for 
specific values of to and tobs corresponding to the Chandra 
observations are given in Table |9] 

Given the observed line widths, it follows from relation 
(l24l) that the nutation broadening can serve as an additional 
5%) correction to the FWHM. 

5. 3. 3 Photoionization and photoexcitation 

So far our analysis was based on the assumption that the 
atoms in the plasma of the jets in SS 433 were ionized and 
excited exclusively by collisions with electrons; in this case, 
the influence of the electromagnetic radiation fields perme- 
ating the jets is negligible. How valid is this assumption? 

As regards the intrinsic X-ray radiation from the jets, it 
should not affect significantly the ionization balance of the 
hot plasma, because the jets are optically thin for the ion- 
izing continuum. However, the gas of the relativistic jets 
in SS 433 can be irradiated by the much more intense X- 
ray radiation produced inside the funnel of a thick accre- 
tion disk near the black hole ( Fabrikal 12004 ). Indeed, if 
SS 433 is similar to the u ltralumin ous X-ray sources ob- 
served in nearby galaxies (lRobertsll2007>) but only turned 
to us by its edge rather than by the accretion disk plane, 
then the X-ray luminosity emitted within the funnels can 
be Lx ~ 10"^^- 10"^^ erg/s (when recalculated to an isotropic 
source). The photoionization fraction of the jet gas in 
the field of such radiation will be determined by the ion- 



ASTRONOMY LETTERS Vol.38 No. 7 2012 



460 



KhabibuUin & Sazonov 



ization parameter ^ = L/nr^ and, depending on the den- 
sity n and the distance to the jet base r, can reach 10^ or 
even greater values. In this case, the influence of pho- 
toionization on the ionization balance in the hot plasma of 
the jets can be significant (e.g., Kallman,McCray ( 1982)X 
Howe ver, Chandra HETGS measurements dMarshall et al.l 



2002b showed that the permitted line dominates in the re- 



solved (in energy) Ka triplets, primarily of SiXIII. This 
means that collisions dominate over photoionization in the 
formation region of these lines. Since the ionization pa- 
rameter ^ = L/nr^ does not change greatly along the X-ray 
jet for both quasi-adiabatic and quasi-cylindrical models, 
collisions dominate over photoionization during the exci- 
tation of levels in the atoms of heavier elements as well. 
Nevertheless, we are going to consider in more detail the 
influence of X-ray radiation from the inner accretion disk 
regions on the ionization and radiative properties of the 
plasma in the jets of SS 433 in a future paper. 

The supercritical accretion disk in SS 433 is also a pow- 
erful source of ultraviolet radiation with a characteristic 
temperature of ^ 5 0000 K and a luminosity of 10^^ erg/ s 
( Dolan et al ] ll997h . This is also confirmed by observations 
of the radio nebula W50 produced by the impact of the jets 
in SS 433 on the interstellar medium ( Fabrikai|2004|). An 
ultraviolet photon of appropriate energy can excite the elec- 
tron from the upper level (^S) of the forbidden transition in 
a helium-like ion t o the upper level CP ) of the intercom- 
bination transition (IPorquet et al.ll200lb . In a strong ultra- 
violet radiation field, this effect leads to a redistribution of 
the radiation fluxes in the forbidden and intercombination 
lines of the triplets in favor of the latter, i.r., to a decrease 
in ratio R, just as in the case of a high gas density. This 
is particularly true f or the lighter elements. Calculations 
(jPorquet et al.ll200l1) (and our Cloudy computations) for a 
blackbody radiation field with a temperature of 50 000 K 
show that the effect for the SiXIII triplet must be significant 
even if the dilution of radiation related to the fact that the 
disk radiation permeates the jets from some distance and at 
some angle is taken into account. Since there is a consider- 
able uncertainty in determining the characteristic temper- 
ature and the geo metry of the ultraviolet radiation source 



( Gies et al.ll 2002h . we will restrict ourselves here only to 



the above remark of a somewhat qualitative character. 

6 COMPARISON WITH OBSERVATIONS 

In this part, we would like to discuss the application 
of our simulation results to analyzing and interpreting the 
Chandra spectra of SS 433 at the p hase of the greates t 
disk opening toward the observer ( Marshall et al.l 120021) . 



The weighted mean line width in these observations cor- 
responded to the scenario of a ballistic flow with an open- 
ing angle 6 = 0.°61 ± 0.°03 ^ 0.01 rad. The jet density 
fie ^ 10 ^"^cm ~^ in the line formation region (7^1 keV) 
was estimated from the SiXIII Ko^ triplet. The physical vol- 
ume of this region, hence, the distance of this region from 



the cone apex (assuming a conical flow with an opening 
angle9 ) was estimated from the emission measure of the 
component with such a temperature in a four-temperature 
model describing satisfactorily the line intensities. Using 
the estimates of the emission measure for the other three 
components and assuming the temperature profile to be adi- 
abatic, the corresponding locations and densities of these 
components were found. As a result, a physical picture 
was obtained that provided satisfactory agreement with the 
observed spectrum but, at the same time, had significant 
shortcomings. Below, we would like to show how this pic- 
ture changes if the scattering effects we considered and 
others (see the "Results" Section) are taken into account. 
Since the triplets of helium-like iron (FeXXV Koc) and sil- 
icon (SiXIII Key) were the main diagnostic tools, we will 
dwell on a comparison of themeasured line characteristics 
with the results of our calculations. 

As has already been noted above, the Chandra HETGS 
spectral resolution (~ 200 @ 6.7 keV) is insufficient for 
investigating the fine structure of the FeXXV Ko, triplet 
(~ 300 @ 6.7 keV is needed). Therefore, when the ob- 
served profile is analyzed, some fixed ratios of the compo- 
nents (f/r and i/r) is assumed (e.g., f/r = 0.23 and i/r = 0.28 
m iNamiki et al.1 (l2003h ). Figures [TTb and [TTh show how 
strongly these ratios are subjected to the scattering effects, 
which points to the necessity of correcting the line widths 
estimated within a simple model. An unresolvable triplet 
with enhanced (compared to the resonance line) forbidden 
and intercombination components can appear as consist- 
ing of broader components with a bright resonance line but 
with a shifted centroid. Bearing in mind the broadening 
by scattering (Fig. [TTb ) and the additional possibilities for 
broadening from the "Results" Section, it should be rec- 
ognized that the models with 6 < 0.01 rad are admissi- 
ble from the viewpoint of observations in the high-energy 
part of the spectrum. In addition, the decrease in triplet 
intensity (Fig. [TTh ) due to the light scattering effects im- 
plies that the emission measure of the hottest parts of the 
jet estimated from the observed iron line intensities and the 
heavy-element abundances estimated from the ratio of the 
observed intensities of the corresponding lines should be 
corrected significantly. Thus, the scattering effects affect 
significantly the jet parameters (6,^eo,^o) estimated from 
the FeXXV Koc triplet characteristics in Chandra observa- 
tions. 

In contrast, the remarkable spectral characteristics of 
the Chandra instruments in the low-energy region allow 
the fine structure of the helium-like silicon triplet (and, to 
some extent, the SXV Ko^ triplet) to be investigated, but, at 
the same time, they also require a more detailed study of 
their formation mechanisms for interpreting the observed 
line characteristics. The observed width of the SiXIII 
triplet components (FWHM ^1900 km/s) corresponds to a 
jet opening angle 6 0.015 radQ. This makes the inter- 



It is worth noting that in lMarshall et al the weighted mean 
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pretation of the widths in terms of quasi-cyUndricalmodels 
very problematic, even if the possible broadening effects 
we considered are taken into account (Fig. [T2b). Thus, 
the models with B > 0.01 rad are much more preferential. 
However, the intensity of the SiXIII Koc triplet decreases 
significantly as a result of scatterings even in this case (Fig. 
[T2h). This implies the necessity of a significant correction 
of the estimated emission measure for the jet region with 
r ~ 1 keV, though slightly less significant than that for the 
hotter region at the jet base. 

Figures [T2b and [T2h present the ratios R and G for the 
silicon triplet derived in our calculations by taking into ac- 
count the scattering effects and the contribution from the 
MgXII Ly^ doublet, whose energy almost coincides with 
the energy of the forbidden silicon line (see also Section 
5.3.1). These results of our calculations are compared 
with the ratios (confidence intervals) R and G measured by 
lMarshalletal.1 (12002): 7?= 1.18±0.26and G = 0.92±0.13. 
The following conclusions can be drawn from this com- 
parison. First, the measured value of G is consistent with 
most of the quasi-adiabatic and quasi-cylindrical jet mod- 
els considered. Second, the measured value of R points to 
ro© ^ 10^ cm, which together with 9 :^ 0.01 rad estimated 
from the line broadening giv es rp ^ 10^^ cm. Th is is consis- 



tent with the estimates from Medvedev. Fabr ika (201 0|) but 



exceeds appreciably the value obtained by Marsha ll et al 
(l2002h - 



^0 



2 X 10 ^^cm, and, at the same time, is less 



than ro ~ 10 ^^cm found by lFilippova et al.l (2006) from the 
eclipse of the jets by the optical companion. 

The presented comparison of the results of our numerical 
calculations w i th the results of observations published by 
Marshall et al.l (|2002|) is more likely qualitative than quan- 
titative in nature. Obviously, comparison of the spectral 
predictions of our model for the scattering of radiation in 
the jets of SS 433 dir ectly with the Chandra HETGS data 
([Marshall etal.ll2002h is required to formulate more reli- 
able conclusions. However, this is beyond the scope of our 
paper. 



(2) The scattering inside the jets, along with the nuta- 
tional motion and the contribution of satellites, lead to a no- 
ticeable additional line broadening. This can lead to over- 
estimates of the opening angle O from the line width in 
Chandra X-ray observations. As a consequence, the mod- 
els with O ~ 0.005 - 0.01 rad remain admissible. There- 
fore, it is worth noting that the more complex scenario of 
a nonco nical relativistic flow (for example, ongoing colli- 
mation; iNamiki et al.l (l2003h ). which can be described by 



the change in 9 along the jet in some approximation, may 
actually be realized. 

(3) There must be broad wings of scattered radiation near 
the brightest lines, distorting the shape of the continuum in 
a certain way. The magnitude of this distortions turns out 
to be of the order of (or even greater than) the contribu- 
tion from the recombination and two-photon components. 
As an illustration, we provide a synthetic broadband spec- 
trum of the approaching jet at a phase corresponding to the 
Chandra X-r ay observations of SS 433 in 1999 at a distance 
of 5 kpc (Mar shaU etaDl2002|) (Fig. [H]). Accurately mea- 
suring the shape of the scattered component can give ad- 
ditional information about the density distribution and the 
temperature profile along the jet. 

(4) The fine structure of the lines is very sensitive to the 
scattering effects (especially in the case of doublets and 
triplets). This makes its investigation a powerful tool for 
diagnosing the jet parameters - the density of the emitting 
region (and, consequently, its size) and the opening angle 
6. New-generation X-ray observatories (primarily Astro- 
H) equipped with spectrometers with a resolution of several 
electronvolts (microcalorimeters) will allow one to quickly 
(in an exposure time of 10 ks, Astro-H will be able to col- 
lect ~600 photons in the FeXXV (6.7 keV) line and, con- 
sequently, ~60 photons in the broad wings) get a clear idea 
of the mechanisms for the formation of X-ray lines in the 
spectrum of the jets in SS 433. 



7 CONCLUSIONS 

Let us briefly list our main conclusions: 

(1) The effect of a decrease in line intensity due to scat- 
tering inside the jet turns out to be very pronounced, but it 
does not exceed 60% in magnitude on the entire grid of pa- 
rameters. Thus, such anomalies as the excess of radiation 
in the region of the helium-like nickel triplet are difficult to 
explain in terms of the scattering effects alone, it is quite 
possible to make the problem less dramatic. At the same 
time, diagnosing the jet plasma parameters based on the 
line intensity ratios turns out to be impossible without al- 
lowance for the scattering effects, because the introduced 
shift can lead to significant systematic errors. 

width FWHM ^1710 km/s measured from all lines corresponds to an 
opening angle B ~ 0.°77 ~ 0.014 rad (which is easy to see from Eqs. (3) 
and (4) of this paper), while the authors provide 6^0. °61 ^0.01 rad. 
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Fig. 10. (a) Influence of the Mg XII Ly^ line on the ratio G = (f + i)/r in 
the Si XIII Ka triplet: the solid and dashed curves are for the cases with 
and without Mg XII Ly^, respectively. The computation was performed 
using ATOMDB v2.0.1. (b) The contribution from the satellites of the 
Fe XXV Ka triplet to the radiation intensity at 6.7 keV: the dashed and 
solid curves are for the cases without and with scattering, respectively; 
red-with satellites, black-without satellites. 
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Fig. 11. Analysis of the simulation results for the FeXXV Ka triplet in the Chandra HEG resolution for the approaching jet with Tq = 20 keV at the 
phase of the greatest disk opening toward the observer. The dashed lines (with the corresponding open symbols) correspond to the unscattered line 
characteristics. 
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Fig. 12. Analysis of our simulation results for the SiXIII triplet in the Chandra HEG resolution for the approaching jet with Tq = 20 keV at the 
phase of the greatest disk opening toward the observer. The dashed lines (with the corresponding open symbols) correspond to the unscattered line 
characteristics. The ratios R and G are presented by taking in to acc ount the contribution from the MgXII Ly^ doublet. The shaded regions correspond 
to R and G measured in Chandra observations (.Marshall etZll2002h . 
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Fig. 13. Broadband simulated synthetic spectrum (solid line) of the approaching jet (z = - 0.078) at a distance of 5 kpc. The continuum radiation 
(dashed line) was computed using ATOMDB v.2.0.1 (the NoLine model). All of the lines from Table^with allowance made for their scattering inside 
the jet but without allowance for the contribution of satellites were included. Magnified fragments of the broadband spectrum in the low-energy (bottom 
left) and high-energy (bottom right) regions. The dash-dotted line indicates the contribution from bremsstrahlung (the brems model in XSPEC vl2.6.0). 
The difference between the bremsstrahlung and NoLine continua is due to the presence of recombination and two-photon components. 
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